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Most of our understanding of bacterial chemotaxis comes from studies of Escherichia coli. How¬ 
ever, recent evidence suggests significant departures from the E. coli paradigm in other bacterial 
species. This variation may stem from different species inhabiting distinct environments and thus 
adapting to specific environmental pressures. In particular, these complex and dynamic environ¬ 
ments may be poorly represented by standard experimental and theoretical models. In this work, 
we study the performance of various chemotactic strategies under a range of stochastic time- and 
space-varying attractant distributions in silico. We describe a novel type of response in which the 
bacterium tumbles more when attractant concentration is increasing, in contrast to the response of 
E. coli, and demonstrate how this response explains the behavior of aerobically-grown Rhodobacter 
sphaeroides. In this “speculator” response, bacteria compare the current attractant concentration 
to the long-term average. By tumbling persistently when the current concentration is higher than 
the average, bacteria maintain their position in regions of high attractant concentration. If the cur¬ 
rent concentration is lower than the average, or is declining, bacteria swim away in search of more 
favorable conditions. When the attractant distribution is spatially complex but slowly-changing, 
this response is as effective as that of E. coli. 

INTRODUCTION 

Movement of Escherichia coli consists of periods of 
running punctuated by tumbling events where the bac¬ 
terium randomly changes direction. This can result in 
successful chemotaxis when the probability of initiating 
a tumble per short time interval (the tumbling rate) is 
a function of the concentration of attractant experienced 
by the bacterium. In the absence of attractant, an E. coli 
bacterium has a constant, basal tumbling rate. When the 
attractant concentration experienced by the bacterium is 
increasing, tumbling events become less frequent, so the 
bacterium has longer runs in the direction of increasing 
attractant. Gonversely, when the bacterium detects a de¬ 
creasing attractant concentration, it tumbles more often, 
shortening its runs down attractant concentration gradi¬ 
ents. In E. coli, this mechanism involves an excitation 
pathway, inhibiting tumbling, and an adaptation path¬ 
way that methylates the receptors, decreasing their sen¬ 
sitivity and thereby attenuating the excitation pathway 
in the continued presence of attractant. This “adaptive” 
response allows the bacterium to locate and stay in re¬ 
gions of high attractant concentration. When attractant 
concentration plateaus, the tumbling rate returns to the 
basal rate, a phenomenon called “perfect adaptation”. A 
consequence of perfect adaptation is that the response is 
independent of the absolute concentration of attractant 
and depends only on differences in concentration experi¬ 
enced by the bacterium. 

While E. coli has been instrumental for our under¬ 


standing of chemotaxis, other bacteria show a consid¬ 
erable variety of chemotactic mechanisms and behav¬ 
iors mu. For example, responses have been identi¬ 
fied which show very little adaptation, such as in certain 
cultures of Rhodobacter sphaeroides [3]- Furthermore, 
aerotaxis in E. coli is thought to involve the Aer recep¬ 
tor which lacks methylation sites |3]) suggesting a lack 
of adaptation. Strangely, some bacteria seem to tum¬ 
ble more in the presence of attractant. This appears 
to be the case in many mutant strains, such as aero¬ 
taxis [S] and redox taxis mutated E. coli, aerotaxis 
in mutated Salmonella typhimurium |5]j s-nd phototaxis 
in mutated Halobacteria [7]- In 5. typhimurium chemo¬ 
taxis, this alternative response can be caused by one of a 
number of single point mutations [5]. Interestingly, this 
behavior was also found in wild-type, aerobically-grown 
R. sphaeroides [S]. This response seems paradoxical, as 
the bacterium might be expected to run in the direction 
of decreasing attractant concentration and tumble more 
when it detects an increase in attractant concentration, 
leading to the accumulation of bacteria away from the 
attractant. In addressing these puzzling results, Gold¬ 
stein and Soyer demonstrated the chemotactic efficacy 
of a non-adaptive “inverted” response [TITl ITT] . With 
this strategy, bacteria respond to the absolute attractant 
concentration, tumbling more and therefore maintaining 
their position in regions of higher attractant concentra¬ 
tion. Although less effective than the adaptive response, 
the inverted response requires only low receptor sensi¬ 
tivity, and could function in the absence of effective re- 
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ceptors by coupling to the cell metabolism m- As we 
show, there are discrepancies between the inverted re¬ 
sponse and the response observed in aerobically-grown 
R. sphaeroides. 

Why do different bacteria exhibit different chemotac- 
tic responses? One possible reason is that different bac¬ 
teria have evolved for different environments. For ex¬ 
ample, while E. coli might be expected to inhabit a 
resource-rich environment, marine bacteria experience a 
harsh environment in which attractant is localized in 
short-lasting patches [12] with attractant concentration 
inside patches being 3 to 6 orders of magnitude higher 
than outside M- This has led to a number of marine- 
specific evolutionary adaptations such as higher running 
speed in Pseudomonas haloplanktis US] and run-and- 
reverse (as opposed to run-and-tumble) chemotaxis in 
over 70% of marine bacterial species m- Unfortunately, 
most experimental and theoretical studies to date con¬ 
sider chemotaxis in response to step functions or simple 
gradients m, providing limited insights to how chemo¬ 
taxis would function in different types of environments. 
In particular, there have few studies [niin] analyzing 
chemotactic strategies in stochastic environments, which 
are likely to be the most important environment during 
the evolution of a chemotactic response. 

In this work we construct a model of a stochastic at¬ 
tractant distribution. By adjusting the manner in which 
attractant concentrations vary in time and space, the at¬ 
tractant distribution can mimic a range of environments 
that one might expect to find in nature. We use this 
model to study how performance and optimal properties 
of various chemotactic strategies vary as a function of 
environmental conditions. In particular, we describe a 
new type of chemotactic strategy called the speculator 
response. It differs from the adaptive response in that 
the tumbling rate increases with increasing attractant 
concentration; furthermore, the bacterium makes tem¬ 
poral comparisons of attractant concentration which dis¬ 
tinguishes the strategy from the inverted response. We 
demonstrate the effectiveness of the speculator strategy 
and its remarkable match to the paradoxical response 
seen in wild-type, aerobically-grown R. sphaeroides |5]. 


RESULTS 

To model chemotaxis, we consider a single bacterium 
in a one-dimensional space with periodic boundary con¬ 
ditions and a distribution of attractant. The bacterium 
can run to the left or right, or tumble, a and /3 denote 
the rates at which the bacterium starts and stops tum¬ 
bling. While P is assumed to be constant, the basal rate 
ao is modulated by the chemotactic response of the bac¬ 
terium to the experienced attractant concentrations. In 


particular, at time t, 

a(t) = max ^0, oo-I- J R(t — dt'^ 

( 1 ) 

where R{t) is the chemotactic response function and 
c(xB(t),t) is the attractant concentration that the bac¬ 
terium experiences at position XB(t) at time t |20j . In 
contrast to m, we do not assume deviations from ao to 
be small. R{t) is represented as {AjT + BtjT"^) exp(—t/r) 
where r controls the memory length, i.e. how far back 
in the past the bacterium “remembers” attractant con¬ 
centrations, and A and B together determine the sensi¬ 
tivity and the characteristics of the response: adaptive, 
inverted, or speculator. In the adaptive response, A and 
B are constrained such that A < 0 and B = —A. This 
gives rise to a response function shown in red in Fig. & 
that has a positive and a negative lobe. The red curve 
in Fig. & illustrates the changes in a due to attractant 
addition and removal when a response function of this 
type is used. When c{xB{t),t) is increasing in time, 
such as when attractant is added, the negative lobe of 
R{t — t')c{xB{t'),t') has a larger area than the positive 
lobe, making the integral in Eq. Q negative, leading to 
an a that is smaller than ao, resulting in a decrease in 
tumbling, as the red curve in Fig.[l]D shows att = 50. The 
opposite happens when attractant is removed, resulting 
in an increase in tumbling at t = 350. The constraint 
B = —A results in equal areas of the positive and nega¬ 
tive lobes of R{t), ensuring perfect adaptation and a basal 
tumbling rate (a = ao) for 50 < t < 350. In the inverted 
response (Figure [^, blue curve), A > 0 and B — 0, lead¬ 
ing to a single-lobe response function. This results in a 
higher tumbling rate in the presence of attractant and a 
lower rate when attractant is absent, as the blue curve 
in Fig. shows. We also investigate a new type of re¬ 
sponse which we name the “speculator” response for rea¬ 
sons explained below. In the speculator response, A > 0 
and B < 0, leading to a double-lobe response function 
(Figure [^, green curve) that looks roughly like the neg¬ 
ative of the adaptive response function (Figure [^, red 
curve). This causes increased (decreased) tumbling when 
c{xB{t),t) is increasing (decreasing) in time (Figure]^, 
green curve). The constraint of perfect adaptation is re¬ 
laxed in the speculator response, so the areas of the pos¬ 
itive and negative lobes are unequal (Figure [^, green 
curve). This allows the steady-state a in the presence 
of constant attractant concentration to be different from 
ao, as shown for times 100 < t < 350 (Figure [^, green 
curve). The double-lobe response functions of adaptive 
and speculator response cause bacteria to make tempo¬ 
ral comparisons of attractant concentration, while the 
lack of perfect adaptation in the inverted and speculator 
responses causes bacteria to respond to absolute attrac¬ 
tant concentrations. Note that Fig. & and Fig. are 
purely illustrative; they do not reflect real or optimized 


3 



t 




0 100 200 300 400 

t 


0 100 200 300 400 

t 


FIG. 1. Illustrative examples of the response function and the corresponding timecourses of a in response to 
attractant step functions. Timecourses of a in response to attractant step functions for optimized responses. 

a) Illustrative examples of the response function for adaptive (red), inverted (blue) and speculator (green) response. The 
parameters used are A = —20, B = 2Q and r = 5 for adaptive, A = 20, B = 0 and r = 5 for inverted, and A = 20, i? = —18 
and r = 5 for speculator response, b) Changes in a in response to step changes in attractant concentration for the responses 
from part a), oq = 10 for all responses. Attractant (concentration of 1) is added at t = 50 and removed at t = 350. c) Changes 
in a in response to step changes in attractant concentration for adaptive (red), inverted (blue) and speculator (green) responses 
optimized for T = 10'* and L = 100 (ao = 0.0084, A = —1500, B = 1500 and r = 0.020 in the adaptive response, olq = 0.0063, 
A = 4.4, B = Q and r = 5.0 in the inverted response, and ao = 0.0089, A = 74, B = —67 and r = 33 in the speculator 
response). Attractant (concentration of 1) is added at t = 50 and removed at t = 350. Note the change of scale on the y-axis. 


responses. 

The one-dimensional virtual world in which the bac¬ 
terium moves contains a stochastic attractant distribu¬ 
tion which varies in both time and space. Two parame¬ 
ters, T (correlation time) and L (correlation length), de¬ 
termine the dynamics of the attractant distribution. T is 
the timescale on which attractant concentrations change, 
while L determines the distance between peaks of attrac¬ 
tant concentration; the shorter L, the more numerous 
and narrow the peaks are and the shorter the distances 
between them. The average amount of attractant avail¬ 
able in the world is independent of T and L. Fig. shows 
how the distribution looks at different combinations of T 
and L. [] illustrates the distribution dynamics as a func¬ 
tion of T and L. 

The framework described above allows us to assess the 
performance of a chemotactic response characterized by 
the response parameters Oq, /?, B and r at a chosen 
combination of attractant distribution parameters T and 
L. Performance, or fitness, of a response, is equal to the 
average cell division rate, which we approximate as the 
inverse of the time it takes the bacterium to experience 
a specified amount of attractant. For any chemotactic 
strategy in any stochastic environment, we can optimize 
the response parameters to maximize the bacterial fit¬ 
ness. Performing this optimization for different strategies 
(by applying appropriate constraints on A and B) under 
different combinations of T and L allows us to explore 
the performance of the different strategies and how this 
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FIG. 2. Examples of the stochastic attractant distri¬ 
bution at different combinations of T and L. Every 
row corresponds to a different combination. The left (right) 
panels show the distribution at time t = 0 (t = 100) in the 
simulation, a) T = 10^, L = 100. b) T = 10'*, L = 20. c) 
T = 100, L = 100. d) T = 100, L = 20. 


performance varies with T and L. 

Figure shows the optimal fitnesses of adaptive, in¬ 
verted and speculator responses as a function of T and L. 
Fitnesses are scaled by the fitness of a non-chemotaxing 
bacterium, whose fitness is independent of T and L: a 
bacterium with a relative fitness of 4 therefore takes 4 
times less time to experience the same amount of attrac¬ 
tant than a non-chemotaxing bacterium. In all strate¬ 
gies, fitness increases with increasing T and decreasing 
L. As T increases, attractant concentrations change more 
slowly, making it easier for bacteria to track attractant 
peaks. At shorter L, fitness is higher because there are 
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more attractant peaks and they are closer to one another. 
This means that if a bacterium loses track of a peak, or 
a peak diminishes in amplitude over time, the bacterium 
only needs to travel a short distance to reach another 
peak. 



FIG. 3. Fitness of optimized adaptive (red), speculator 
(green) and inverted (blue) response as a function of 

T and L. The two panels show different views of the same 
plot. For each chemotactic strategy and combination of T 
and L, fitnesses are averaged over the last 600 generations of 
up to 3 replicate simulations. 

In addition to the previously studied adaptive and in¬ 
verted responses, we characterize a novel chemotactic 
strategy. This “speculator” response, despite its seem¬ 
ingly paradoxical nature, is more fit under all studied 
conditions than the inverted response, although less fit 
than the adaptive response. Interestingly, at T = 10^, 
L = 20, the fitness of the speculator response is nearly 
identical to the fitness of the adaptive response. To un¬ 
derstand the mechanism of the speculator response, we 
consider the optimal values of response parameters (Ta¬ 
ble |T]). The lack of perfect adaptation (optimal |A| > \B\) 
means that the bacterium will more often start to tum¬ 
ble when the attractant concentration is high, as shown 
in green in Fig. [T]:; the low value of /3 means that the 
bacterium will then continue tumbling, remaining in the 
region of high attractant. Consequently, the specula¬ 
tor response, like the inverted response (Figure [^, blue 
curve), results in frequent long tumbles at high attrac¬ 
tant concentrations. In contrast to the inverted response, 
the double-lobe response function of the speculator re¬ 
sponse results in a tumbling rate sensitive to the rate of 
change of the attractant. The long memory of the spec¬ 
ulator response (large r) allows sensitivity to long-term 
trends; this sensitivity, combined with the double-lobe 
response function, results in two important dynamical 
properties. Firstly, the bacterium compares recent at¬ 
tractant concentrations with a long-term average, tum¬ 
bling more when the recent past is more favorable than 
the average, and therefore maintaining its position in re¬ 
gions of higher attractant concentration. Secondly, the 
bacterium is able to sense improving and worsening con¬ 
ditions at its current location. In particular, a decline in 
the attractant concentration results in a decrease in a. 


allowing the bacterium to swim away from a peak when 
conditions are changing for the worse. Swimming away 
leads to a further decrease in a, setting a feedback loop 
in motion, resulting in continued swimming until a new 
optimum is reached. The speculator response is therefore 
analogous to the behavior of investors in financial mar¬ 
kets: when the current performance is lower than the 
average, or when investment values are falling, specula¬ 
tors seek higher returns by abandoning their current po¬ 
sition and investing elsewhere—hence the name “specula¬ 
tor” response. The behavior of the speculator response, 
compared with the adaptive and inverted responses, is 
illustrated infl 

Significantly, the time course of a in the speculator 
response closely matches the time course of probabil¬ 
ity of tumbling in aerobically-grown R. sphaeroides (see 
Fig. 2A in i). Fig.|^shows a curve-fit of our model of the 
speculator response to the digitized data of [5]. The close¬ 
ness of the fit provides strong evidence that aerobically- 
grown R. sphaeroides uses the speculator response to re¬ 
spond to Na-succinate. Experimental results show that 
aerobically-grown R. sphaeroides performs well in swarm 
plates [S], demonstrating the efficacy of this response. 



Time (s) 


FIG. 4. Curve-fit of our model of the speculator 
response (green curve) to the digitized data of [9] 
(black circles). The response is described by the follow¬ 
ing response parameters: oo = 0.074s“^, (3 = 0.034s“^, 
A = 1300mM“^s“^, B = —1000mM“^s“^ and r = 71 s. At¬ 
tractant concentration is set to 0.001 mM in correspondence 
to |5]. The curve-fit is obtained by optimizing the response 
parameters of the speculator response to minimize the least- 
squares fit between the model and the digitized data. The 
data are digitized in MATLAB using the function imfindcir- 
cles |21) . 

As Fig. shows, at T = 10^, L = 20, fitnesses of 
the adaptive and speculator responses are very similar 
despite the different mechanisms behind their chemotac¬ 
tic strategies. To better understand these strategies, we 
create a simple attractant distribution which consists of 
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TABLE I. Range of optimal response parameters at T = 10^, L = 50. 3 replicate simulations are run for each 

chemotactic strategy at T = 10^, L = 50. For a given strategy and response parameter, mean values of the parameter are 
calculated separately in each replicate simulation by averaging the values of the parameter over the last 600 generations. Each 
range in the table is composed of the lowest and highest mean values obtained. 



T 

oo 

P 

A 

BjA 

Adaptive 

(0.072, 0.27) 

(0.012, 0.062) 

(10, 47) 

(-11000, -5700) 

-1 

Inverted 

(5.6, 8.3) 

(0.0024, 0.0049) 

(0.098, 0.68) 

(2.0, 6.2) 

0 

Speculator 

(39, 45) 

(0.0030, 0.046) 

(0.086, 0.44) 

(47, 96) 

(-0.86, -0.82) 


two Gaussians (at positions 25 and 75 in a world with a 
length of 100) oscillating in amplitude out of phase with 
each other: when one Gaussian is at full amplitude, the 
other has amplitude of zero. Amplitude, period of oscil¬ 
lation and width of the Gaussians are roughly matched 
to T = 10'*, L = 20 of the stochastic attractant distri¬ 
bution. For each of the two chemotactic strategies, we 
take a bacterium optimized for T = 10*, L = 20 of the 
stochastic attractant distribution and simulate its move¬ 
ment in the virtual world with the two Gaussians. Fig. 
shows the mean position of the bacteria as a function 
of 9, the phase of the oscillations. Between 9 = 0 and 
9 = 1250 (the period is 5000) the left Gaussian at posi¬ 
tion 25 is higher than the right Gaussian at position 75, 
but is decreasing. In the adaptive response (red curve), 
the bacterium is close to the top of this Gaussian dur¬ 
ing this period. The bacterium shows little movement 
toward the right Gaussian at position 75 until the right 
Gaussian is significantly higher i.e. 9 > 1250. In the 
speculator response (green curve), the bacterium cannot 
track the top of the Gaussian as well as in the adaptive 
response, as evidenced by the large standard deviation 
around position 25 (green shading). However, the bac¬ 
terium more quickly adapts to the changing attractant 
levels, leaving the declining left Gaussian and moving to¬ 
wards the growing right Gaussian sooner. 

The strengths of the adaptive and speculator responses 
therefore lie in exploitation and exploration, respectively. 
In the adaptive response, the bacterium can track the 
top of a peak efficiently while in the speculator response, 
the bacterium is better at leaving the declining peak and 
finding the increasing peak. The exploitation behavior of 
the adaptive response is analogous to a hill-climbing algo¬ 
rithm, which efficiently finds, but may get stuck at, a lo¬ 
cal optimum. The exploration behavior of the speculator 
response is more analogous to a Monte Garlo search algo¬ 
rithm in that the bacterium may leave a peak in search 
of a higher peak at the cost of its ability to track the 
peak top efficiently. This explains the trend in Fig. 
for large L, the number of attractant peaks is small, and 
exploiting a given peak is more important than explor¬ 
ing new peaks. Under these conditions, the adaptive re¬ 
sponse is significantly more effective than the speculator 
response. At short L, there are multiple peaks in the 
environment, each of which has a different amplitude. 
Under such conditions, the exploration behavior of the 



FIG. 5. Mean position of bacteria performing adap¬ 
tive (red curve) or speculator (green curve) response 
as a function of 9, the phase of the oscillations. Shad¬ 
ing shows the standard deviation of the position. The Gaus¬ 
sians are centered at positions 25 (dashed gold line) and 75 
(dashed turquoise line) and have a standard deviation of 3. 
Amplitudes of the Gaussians (not to scale, maximum ampli¬ 
tude is 1) are shown as a function of 9 in the bottom part 
of the figure for the Gaussian at position 25 (gold curve) and 
75 (turquoise curve). The parameters used are oq = 0.0084, 
P = 54, A = —1526 and r = 0.020 for the adaptive and 
no = 0.0089, P = 0.056, A = 74, R = -67 and r = 33 for the 
speculator response. 


speculator response allows the bacterium to locate higher 
peaks, while the exploitation behavior of the adaptive re¬ 
sponse may lead to the bacterium tracking a suboptimal 
peak. At T = 10"*, L = 20, the two strategies are ap¬ 
proximately equally effective, giving rise to very similar 
fitnesses (Figure [^. 

We next consider the optimal values of the response 
parameters. In the adaptive response, r (the memory 
length) is very short (Table , allowing the bacterium 
to quickly adjust to small displacements from attractant 
optima. P, the rate at which the bacterium stops tum¬ 
bling, is quite large, corresponding to short-lasting tum¬ 
bles characteristic of chemotaxis in E. coli [^. High sen¬ 
sitivity (large |A| and |H|) is necessary for the bacterium 
to respond to small differences in attractant concentra¬ 
tion characteristic of small displacements from the top 
of an attractant peak. High sensitivity is responsible for 
the high a when the attractant is removed at t = 350 
in Fig. (red curve). Optimal uq, the tumbling rate in 
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the absence of attractant (or under constant attractant 
in case of perfect adaptation), is very low in all strate¬ 
gies, as is evident from Fig. [T]:. Low oq enables bacteria 
to run persistently in order to find regions with more fa¬ 
vorable conditions more quickly. The near-zero value of 
Qfo removes the possibility of a going below oq, eliminat¬ 
ing the response to increasing attractant in the adaptive 
response (Figure [^, red curve, t = 50). 

In the inverted response, the bacterium tumbles more 
at higher concentrations of attractant (Figure [^, blue 
curve). T is longer than for the adaptive response, al¬ 
lowing the bacterium to integrate over short-term fluc¬ 
tuations. In both inverted and speculator responses, /3 
is much lower than in the adaptive response, resulting in 
significantly longer tumbles. This is central to the strate¬ 
gies, as it is the persistence of position when tumbling 
that allows bacteria to stay in regions of high attractant 
concentration. The sensitivity is lower than in the adap¬ 
tive response, in agreement with simple models showing 
that the inverted response is optimized by lower sensitiv¬ 
ity m- Sensitivity needs to be tailored to the range of 
attractant concentrations the bacterium experiences: if 
it is too low, the bacterium will run past high concentra¬ 
tions of attractant; if it is too high, the bacterium will 
tumble at low concentrations of attractant, never reach¬ 
ing higher concentrations. 

DISCUSSION 

In this work we describe a new chemotactic strategy, 
termed the speculator response, in which the bacterium 
compares the current attractant concentration with a 
long-term average; if the current concentration is higher 
than the long-term average, the bacterium tumbles per¬ 
sistently to maintain its position. On the other hand, 
declines in the current concentration will increase the 
probability that the bacterium will swim away to a higher 
peak. By considering stochastic attractant distributions, 
we show that under slowly-changing but spatially com¬ 
plex attractant concentrations (large T, small L), the 
speculator response is almost as efficient at co-localizing 
with attractant as the adaptive response of E. coli (Fig¬ 
ure 1^. While the adaptive response achieves high fitness 
by accurately tracking the top of an attractant peak, the 
speculator response enables the bacterium to explore its 
environment and find higher peaks more efficiently (Fig¬ 
ure]^. 

The speculator response closely matches the response 
observed in wild-type, aerobically-grown R. sphaeroides 
(Figure]^. The optimized response parameters from our 
simulations are in arbitrary units, and cannot be directly 
compared with those obtained by the fit to the wild-type 
response (Figure]^. Interestingly, however, the ratio of 
B to A which quantifies the extent of departure from 
perfect adaptation {BjA = —1 corresponds to perfect 


adaptation) is similar between the optimized values ob¬ 
tained from the simulations and the response observed 
in aerobically-grown R. sphaeroides (—0.86 and —0.82, 
respectively). Furthermore, we can acquire a rough esti¬ 
mate of the ratio of ts /ta (where ts and ta are the values 
of T in the speculator and adaptive responses) by com¬ 
paring the values of r in aerobically-grown R. sphaeroides 
(Figure]^ and wild-type E. coli [IS]- This ratio (71/1 = 
71) is of similar order of magnitude to the ratio for the op¬ 
timized simulated responses (43/0.20 = 215), despite the 
multitude of differences between wild-type R. sphaeroides 
and E. coli. 

The optimized adaptive response possesses high sensi¬ 
tivity (large |A| and \B\; Tablejl]) consistent with experi¬ 
mental results from E. coli |2Tj. Furthermore, (3, the rate 
at which the bacterium stops tumbling, is high, which 
is in line with the short tumbles observed in real bacte¬ 
ria m- In contrast to real bacteria, the optimized bacte¬ 
ria have a lower ao, and thus tumble less than real bac¬ 
teria when attractant concentration is increasing (Fig¬ 
ure [^, red curve, t = 50). This may be an artifact of 
modeling chemotaxis in a one-dimensional environment: 
in a three-dimensional environment, tumbling may assist 
the bacterium in finding even steeper paths to attractant 
optima. 

Our model does not take into account the motility- 
associated energy costs of the different chemotactic 
strategies. For instance, R. sphaeroides does not actively 
tumble, but rather stops running and lets rotational dif¬ 
fusion generate the re-orientation, reducing the costs of 
strategies that involve longer tumbles m- The specu¬ 
lator response therefore might have emerged partly be¬ 
cause R. sphaeroides uses rotational diffusion to achieve 
tumbling. Alternatively, rotational diffusion might have 
emerged in response to the bacterium using a strategy 
that involves long tumbles. In addition, by necessity we 
are confined to a relatively small range of T and L; other 
conditions might exist (such as larger T and smaller L) 
that would favor the speculator response even more. 

Our approach differs from that of other studies in 
that we consider realistic attractant distributions and 
extended tumbling times. The latter is essential for the 
speculator response to work as it allows bacteria to main¬ 
tain their position in regions of high attractant concen¬ 
tration. Previous studies nnnunzi modeled tumbles as 
instantaneous after chemotaxis in E. coli m, however, 
experimental evidence from other bacterial species shows 
longer tumbling times [251 US] • Our results add to the 
growing body of evidence that extended tumbles allow 
for emergence of other modes of chemotaxis [iniiniiss. 

Most studies to date considered chemotaxis in response 
to step functions or simple gradients [T7]. While this 
is important for our understanding of the basic mech¬ 
anisms of chemotaxis, we should recognize that chemo¬ 
tactic strategies were inevitably shaped by the environ¬ 
ments the bacteria inhabited. For example, studies in 
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marine bacteria unearthed specific adaptations to ma¬ 
rine environments US na ng, highlighting the need to 
study chemotaxis in the context of realistic attractant 
distributions. Here, we propose a model of a stochastic 
attractant distribution which allows us to compare the 
performance of various chemotactic strategies under dif¬ 
ferent environments and study how optimal properties 
of chemotactic responses change as a function of envi¬ 
ronmental conditions. This can also help us characterize 
the environmental conditions based on the strategies that 
have evolved. Further characterization of natural envi¬ 
ronments m will allow theorists to construct more de¬ 
tailed attractant distributions and advances in microflu¬ 
idics technologies will enable these environments to be 
reconstructed in laboratory settings m- 


METHODS 

Stochastic attractant distribution 


We generate our stochastic attractant distribution by 
summing over cosine and sine modes with different mode 
numbers p so that the concentration at position x and 
time t along the virtual world is calculated as 


c{x,t) = max 0, y: Xp{t) cos Cp + Yp{t) sin ^p (2) 


p=i 


where Xp{t) and Yp{t) are stochastic weights, I is the 
length of the one-dimensional virtual world {I = 100), 
^p = 2'iTpx/l and p* = l/L is the largest mode included 
in the sum above. Wp(0) = Fp(0) = 0 for all ps, and are 
updated at intervals of Ate = T/lOO according to: 


Xp{t + Ate) — ^p(f)(l ~ Ate/T) + rjp{t) 



( 3 ) 


where ? 7 p(t) is a white noise Gaussian random process 
= ^it~t')^pq)j generated by a random num¬ 
ber sampled from a normal distribution with mean 0 and 
variance 1. A similar expression is used for Yp{t + Ate). 
By construction, this results in a Markov process with 
correlation time T and approximate correlation length 
L = IIp*. 


Chemotaxis 

The attractant distribution is equilibrated for a period 
of at least T. Before a bacterium is introduced, a is ini¬ 
tialized based on the equilibrated attractant distribution. 
The bacterium is then released and the state of the bac¬ 
terium (whether it is running or tumbling) is updated ev¬ 
ery Ats = min(r, L/v, t/ 20) where v is the speed of the 


bacterium when running (r; = 1). A Monte Carlo scheme 
is used to decide whether the bacterium starts tumbling 
(running) given that it was running (tumbling) previ¬ 
ously, assuming first-order dynamics of a 2-state system. 
When the bacterium stops tumbling, it starts running 
left or right with equal probability. 


Optimization 


Mutagenesis followed by selection constitute one gen¬ 
eration of the optimization. In the first generation, all 
response parameters are initialized randomly from a uni¬ 
form distribution between 0 and 1 (but see below). In 
the adaptive response, only B is initialized and mutated, 
A is set to —B (at T = 10"*, B is initialized between 
1 and 10). In the speculator response, B is initialized 
randomly between 0 and —1. In every generation, one 
response parameter is chosen at random and mutated. 
Parameters are mutated on a log scale by a transforma¬ 
tion exp(logg(a) -|-r) where a is the parameter being mu¬ 
tated and r is a random number sampled from a uniform 
distribution between —0.2 and 0.2. Further constraints 
on response parameter values are imposed for reasons of 
computational tractability: ag > 10“^ in adaptive and 
inverted response, A > exp(—1) and \B\ > exp(—1) in 
speculator response, r > 0.01 in adaptive response. 

After mutagenesis, the fitnesses of responses described 
by the wild-type and mutant response parameters are de¬ 
termined. This is achieved by letting 10 identical wild- 
type and 10 identical mutant bacteria explore the virtual 
world with the stochastic attractant distribution. Each 
of the 10 wild-type bacteria is subjected to an attractant 
distribution initialized with a different random seed; the 
attractant distributions are then re-used for the 10 mu¬ 
tants. (As the attractant distribution is stochastic, esti¬ 
mates of response fitness are stochastic too. This scheme 
of competing the wild-type and mutant with the same 
attractant distributions is thus used to ensure that lucky 
mutants do not fix.) Each of the bacteria is run until 
it experiences 50T attractant units. (7?m,i) denotes 
the time it took the f-th wild-type (mutant) to experience 
the specified amount of attractant. Fitness of response 
k = {w,m}, Efe, is then calculated as T{j^)i averaged 
over the bacteria. 

Once fitnesses are determined, the probability of ac¬ 
ceptance of the mutation, p^, is calculated using the 
Metropolis-Hastings algorithm: = 1 if and 

Pm = exp((F’m — Fw)/(t/Fw)) otherwise. U, the temper¬ 
ature, is constant at 0.005. Simulations are run until the 
fitness stops increasing and stays constant for at least 
600 generations. 3 replicate simulations are run for each 
chemotactic strategy and combination of T and L. 




SUPPORTING INFORMATION 

SI Video 

Dynamics of the attractant distribution for a) 
T = 10000, L = 100, b) T = 10000, L = 20, c) T = 100, 
L = 100, a) r = 100, L = 20. 


S2 Video 

Dynamics of optimized adaptive (a), inverted 
(b) and speculator (c) strategies under T = 10000 
and L = 20. The parameters used are: Oq = 0.0065, 
P = 4.2, A = —2100, B = 2100, r = 0.016 in adaptive, 
ao = 0.0016, P = 0.048, A = 3.7, B = 0, r = 5.8 in 
inverted and ao = 0.0089, P = 0.056, A = lA, B = —67, 
T = 33 in speculator response. In the adaptive response, 
the bacterium swims up attractant gradients and tum¬ 
bles when it experiences a decrease in attractant con¬ 
centration. This leads to an oscillatory behavior around 
peak maxima. In the inverted response, tumbling rate 
increases with increasing attractant concentration. Re¬ 
sponse sensitivity is optimized such that the bacterium 
tumbles most persistently at attractant concentrations 
which correspond to typical values at attractant maxima. 
However, this means that the bacterium can get stuck at 
sub-optimal concentrations on large peaks. The specu¬ 
lator response compares the current concentration of at¬ 
tractant with a long term average. If the current concen¬ 
tration is greater than this average, the bacterium tum¬ 
bles more. If the current concentration is lower than the 
average, or declining, the bacterium swims away, leaving 
the peak to search for higher attractant concentrations. 
The bacterium will typically run past peaks if their am¬ 
plitude is lower than the peak it just left. 
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